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FINITE ELEMENT APPROXIMATION OF THE PARABOLIC 
FRACTIONAL OBSTACLE PROBLEM* 

ENRIQUE OTAROLAt AND ABNER J. SALGADO* 

Abstract. We study a discretization technique for the parabolic fractional obstacle problem 
in bounded domains. The fractional Laplacian is realized as the Dirichlet-to-Neumann map for a 
nonuniformly elliptic equation posed on a semi-infinite cylinder, which recasts our problem as a 
quasi-stationary elliptic variational inequality with a dynamic boundary condition. The rapid decay 
of the solution suggests a truncation that is suitable for numerical approximation. We discretize the 
truncation with a backward Euler scheme in time and, for space, we use first-degree tensor product 
finite elements. We present an error analysis based on different smoothness assumptions. 

Key words, obstacle problem, thin obstacles, free boundaries, finite elements, fractional diffu¬ 
sion, anisotropic elements. 
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1. Introduction. In this work we shall be interested in the design and analysis 
of a finite element approximation of the so-called parabolic fractional obstacle problem. 
Let n be an open and bounded subset of with d>l. Given s G (0,1), an obstacle 
-0 : n —)■ R, an initial datum Uq : R and a forcing term f : 17 x (0, T] —>■ R, 

the parabolic fractional obstacle problem asks for a function u : 17 x [0, T] —>■ R that 
satisfies the evolution variational inequality: 

min { dtu -I- (—A)'*u — f, u — ^} = 0 (1.1) 

and U|t=o = Uq- Here, (—A)® denotes the fractional powers of the Laplace opera¬ 
tor, supplemented with homogeneous Dirichlet boundary conditions, which for con¬ 
venience we will simply call the fractional Laplacian. We refer to go for a precise 
definition. We must immediately remark that although our exposition is for the 
fractional Laplacian, our techniques and results are equally applicable to fractional 
powers of a symmetric and uniformly elliptic second order differential operator L, sup¬ 
plemented with homogeneous Dirichlet boundary conditions: Lw = —V-(AVw) +cw, 
with A G GL(R'^)) symmetric and positive definite and 0 < c € C'°’^(i7,R). 

The only caveat is that, at the time of this writing, no regularity results are available 
for (11.11) with the fractional Laplacian replaced by L®. 

The study of numerical techniques for nonlocal problems is a rapidly growing field 
of research. Fractional diffusion has received a great deal of attention in diverse areas 
of science and engineering such as mechanics [3] , biophysics [6] , turbulence [9] , image 
processing [14], peridynamics [15] and nonlocal electrostatics m- In particular, the 
study of constrained minimization problems such as the parabolic fractional obstacle 
problem (ED has received considerable attention. This type of problems arises, for 
instance, in financial mathematics as a pricing model for American options. The 
function u represents the rational price of a non perpetual American option where 
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the assets’ prices are modeled by a Levy process, and the payoff function is ip] see 

[HIMIISI]. 

Based on the Caffarelli-Silvestre extension [5] in previous work we provided a 
comprehensive analysis of the discretization of the linear elliptic case |26) . evolution 
equations with fractional diffusion and Caputo fractional time derivative [25) and the 
elliptic fractional obstacle problem [24] . In this work we proceed in our research pro¬ 
gram and show the flexibility of the ideas developed in |26j by studying the parabolic 
fractional obstacle problem (HH). To the best of our knowledge, this is the first work 
that addresses the numerical approximation of this problem. 

Our presentation is organized as follows. The notation and functional setting 
is described in Section [51 where we also briefly describe, in 112.11 the definition of 
the fractional Laplacian, its localization via the Caffarelli-Silvestre extension ( 112.21) 
and the well-posedness of the fractional parabolic obstacle problem; see 112.31 The 
numerical analysis of problem begins in Section [31 where we discuss a domain 
truncation that allows us, in subsequent sections, to consider a space discretization 
using first-degree tensor product finite elements. The time discretization and its error 
analysis is described in Section (4] The space discretization and its analysis is the 
content of Section [S] we provide an error analysis with minimal ( 1I5.2|) and maximal 
fl l5.3|) regularities. This analysis relies on the construction and approximation prop¬ 
erties of a positivity preserving interpolant. For s > 3/8 we construct an interpolant 
with the requisite properties in Section [ 6 ] 

2. Notation and preliminaries. In this work If is a convex bounded and open 
subset of R'* (d > 1) with polyhedral boundary. Our ideas are equally applicable 
to domains with curved boundaries, but the exposition becomes rather cumbersome 
and so we prefer to avoid it. We will follow the notation of [26] and define the semi- 
infinite cylinder and its lateral boundary by C = O x (0,oo), BlC = dV, x (0, oo). 
For fT > 0 we define the truncated cylinder Cy = O x (0, IT) and its lateral boundary 
dLCy = dQ X (0,fk). We also define the Dirichlet boundary F£> = dLCy U fl x {‘X}. 
Since we will be dealing with objects defined on R'’* and R‘’*+^, it will be convenient 
to distinguish the d -|- 1-dimension. For x G we denote 

X = (x^, ■ ■ ■ ,x^, = (x', = (x', y), x' G R'^, y G R. 

Whenever is a normed space we denote by || • ||x its norm and by X' its dual. 
For normed spaces X and Y we write X ^ Y to indicate continuous embedding. 
We will follow standard notation for function spaces [U [35]. In addition, for an 
open set D C R^, iV > 1, if w is a weight and p G (l,cx3) we denote the Lebesgue 
space of p-integrable functions with respect to the measure uidx by LP{uj,D)\ see 
[ISlEIlEs]. Similar notation will be used for weighted Sobolev spaces. If T > 0 and 
(/ : £> X [0, T] ^ R, we consider as a function of t with values in a Banach space X , 
p : [0,T] 3 1 1 -^ p{t) = p{-, t) G X. For 1 < p < oo we will say that p G T^(0, T-,X) if 
the mapping t >->■ ||^i(t)||x is in LP{0,T). We also introduce the space BV{0,T;X) of 
X-valued functions of bounded variation jS] Definition A.2] 

Yaixg ■■= sup \ V \\g{rj) - g{rj-i)\\x > < oo, 

^ li=i J 

where the supremum is taken over all partitions V = {0 = ro < ■ ■ ■ < rj < ... < rj = 
T} of the time interval [0,T]. We recall that if 5 G BV{0^T; X), then at every point 
to G [0,T) there exists the right limit g^{to) = limtj^tj, g(t) 0 Lemma A.lj. 
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Let /C € N be the number of time steps. We define the uniform time step as 
T = T/IC and we set tfe = fcr, k = 0,. .. ,/C. Given a function w : [0,T] —>■ X, we 
denote = w{tk) G X and w'^ = C X. For any sequence C X, we 

define the piecewise constant interpolant W'^ G L°°{0,T; X) by 

W^t) = w'^+^ tG ( 4 , 4 + 1 ], k = 0 ,...,X-l. 

We also define the piecewise linear interpolant W'^ G C{\p,T]]X) by 

W"(t) = tG [ 4 , 4 + 1 ], fc = 0 ,...,/C-l. 

T T 

The first order backward difference operator c) is defined by . 

We note that dtW'^{t) = for all t G (4-, ^fc+i) and fc = 0,..., /C —1. Finally, 

we also notice that, for any sequence W'^ C X and p G [1, 00 ) we have 



and ||kF^||Loo(o,T;X) = IIW^'^||l~(o,T;Js:) = max{||VF'=||^ : fc = 0 ,...,/C}. 

The relation a < 6 indicates that a <Cb for a constant that does not depend on 
either a or b, but it might depend on the problem data. The value of C might change 
at each occurrence. 

2.1. The fractional Laplacian. For a bounded domain there are several ways, 
not necessarily equivalent, to define the fractional Laplacian; see |26j for a discussion. 
As in [26] we will adopt that based on spectral theory |4]. Namely, since — A : T>(A) C 

—>■ is an unbounded, positive and closed operator with dense domain 

21(A) = 22o(n) n and its inverse is compact, there is a countable collection 

of eigenpairs {A;, (pi}i^n C R"*" x Hq{V,) such that is an orthonormal basis of 

L'^{n) and an orthogonal basis of If re G C^{n) 

w = '^wi(pi, wi= w(fidx', 
ieN 

then, for any s G (0,1), we define {—AYw = is well known, 

the theory of Hilbert scales presented in [231 Chapter 1] shows that V ((—A)®/^) = 
]HI®(H) = [L^(H), i2g (fl)]s, i.e., the real interpolation between L^{n) and i2g (H). Con¬ 
sequently, the definition of (—A)® can be extended by density to the space IHI®(H). If, 
for 0 < s < 1, we denote by IHI“®(H) the dual space of IHI®(H), then (—A)® : IHI®(fl) — 
]HI“®(H) is an isomorphism. 

2.2. The Caffarelli-Silvestre extension problem. The Caffarelli-Silvestre 
result [H1I2SI, requires us to deal with a nonuniformly elliptic equation. With this in 
mind, we define the weighted Sobolev space 

-^i(2/“,C) = {w G fc2^(y“,C) : w = 0 on a+C} . 

Since a G (—1,1), \y\°' belongs to the Muckenhoupt class A 2 (R"'''‘^); see [T2l|33]. Con¬ 
sequently, is a Hilbert space, and smooth functions are dense in Yl\{y°',C) 

(cf. |3S| Proposition 2.1.2, Corollary 2.1.6]). 
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As |26[ (2.21)] shows, the following weighted Poincare inequality holds: 

lkllL^(y-»,c) < l|V«|U2(,.,c), yweHl{y‘^,C). (2.1) 

Then, the seminorm of W{y°‘,C) is equivalent to the norm in H\{y°‘,C). For w G 
H^{y°‘, C) we denote by tro w the trace of w onto fl x {0}. We recall ([^ Prop. 2.5]) 

trnHl{y°‘,C) =M%n), || trn w||H-(n) < (2.2) 

The seminal work of Caffarelli and Silvestre showed that the operator 

(—A)® can be realized as the Dirichlet-to-Neumann map for a nonuniformly elliptic 
boundary value problem. Namely, if G Hj^{y°‘,C) solves 

- V- (y“VW) = 0 in C, = 0 on BlC, d'^U = d^/ on fl x {0}, (2.3) 

where a = 1 — 2s, d'^U = — lim.y_i,o y^^Uy and ds = 2“r(l — s)/r(s) is a normalization 
constant, then u = tiaU G IHI®(r2) solves 

i-AYu = f. (2.4) 

The reader is referred to [26] for a detailed and thorough exposition on how the 
groundbreaking identity given by (12.3|) can be used to design and analyze an efficient 
finite element approximation of solutions to dill). 

2.3. The parabolic fractional obstacle problem. Given an obstacle ip that 
satisfies ip G ^(n) n C{Cl) and -0 < 0 on 917, let 

A1(17) = {w g IHI®(f7) : w{x') > ip{x') a.e. x' G 17} 

be the convex set of admissible functions, and let Ind 5 {-(f 2 ) be its indicator function, 
which, since 1^(17) is closed and convex, is a non-smooth lower-semi-continuous convex 
function. We consider the energy 

+ lnd:^(n) (<(>)• 

With this notation, the parabolic fractional obstacle problem can be understood as 
the gradient flow for J, or an evolution equation for a maximal monotone operator: 
Given an initial datum Uq G L^(17) and a function f G L^(0,T;i^(17)), find u such 
that u(0) = Uq and it solves the differential inclusion 

dtu(7)-I-9J(u(7)) 9 f(7) a.e. 7 g(0,T). (2.5) 

Problem 1)2.51) can be equivalently understood as an evolution variational inequal¬ 
ity: Find u such that u(t) G (^(17) and for a.e. t G (0,r) and all (p € (^(f^) 

(dtu(7),u(7) - (^)z, 2 (o) -p ((-A)ffi(7),u(7) - (p) < (f(i),u(i) - p)L^n), (2.6) 

and u(0) = Uq. Here and in what follows (•, •)L2(n) denotes the inner product of T^(17) 
and (•, •) the duality pairing between IHI“'’(f7) and ]HI'*(f7). 

From these formulations existence and uniqueness of solutions and a priori esti¬ 
mates can be easily obtained with standard techniques on maximal monotone opera¬ 
tors [Sj. For instance, if u solves (12.6L then it satisfies the energy estimate 

l|u|li“(0,T;L2(n)) + l|u|lL2(0.T;H*(n)) ~ 
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where we denoted 

® =2) (uq, f, V”) = l|uo|lL2(n) + IKIlL2(o,T;L2(n)) + ||'*A|Ihs(0)- (2-7) 

If f € L^{0,T; L^(fJ)) and uq € then there exists a unique strong solution, that 

is u G C°([0, T]; L^(r2)) which is locally absolutely continuous in (0,T) and satisfies 
(12.61) at almost every point. In addition, we have that u G H^{0,T; L'^{il)) and the 
mapping t i—>■ J(u(t)) is locally absolutely continuous in (0,T], which implies that 
u G and that the following estimate holds [SJ Theorem 3.6] 

I|dtu(t)||i 2 (n) + dt||u(t)||2,(j.j^ = (f, dtu(t))i 2 (n) a.e. t G (0,r). 

If, moreover, Uq G 'D{dJ) and f G BV{0,T; then u G T]; T^(n)) O 

Proposition 3.3]. Finally, if f G IF^’^(0, T; L^(n)) we have 

idt||dtu(t)||i2(n) < (ft, dtu(t))i2(f2), 
in the distributional sense. 

Let us now use the Caffarelli-Silvestre extension detailed in il2.II to write an 
obstacle problem that is equivalent to (1^ . To do this, we define the set 

tKiC) = G H\{y°‘,C) : trnui(a:^) > ip{x') a.e. x G n| . 

Problem (12.61) can then be equivalently stated as: Find : [0,T] “KiC) such that 
for a.e. t G (0,r) and every (j) G tK{C) 

(trn dt'^(t),tro('^(t) - 4>))L^(a) + a{^{t),^{t) - (j)) 

< (f,trn('^(t) - (^))l 2 (o), (2 .8) 

with tro (0) = Uq. Here the bilinear form a is defined by 

a(w,(()) = ^ y 2 /“VwV(/)da;'dy, Vw,(j) € Hl{y°‘,C). 

The description of the functional setting for problem (12.8|) . together with existence 
and uniqueness results, follow the analysis developed for problem (12.5|) : see [5l[28]. In 
particular, we have the energy inequality 

\\t^n^\\l<^(0,T-,L^{n)) + ll'^llL2(0.T;#i(y°,C)) - 

where 2) is defined in (1^ . We shall be more specific on the smoothness of the 
data and the consequences on the regularity of the solution when we perform the 
discretization and its analysis. Let us now contempt ourselves with mentioning that, 
provided ^ is sufficiently smooth, the following complementarity system holds: 

Z := + ds tro dt^ - d^f >0, tro - '0 > 0, Z (trn - tp) = 0. 

3. Truncation. The variational inequality (|2.8I) is posed on a infinite domain 
and, consequently, it cannot be directly approximated with finite element-like tech¬ 
niques. A first step towards the discretization is to truncate the domain C to a 
bounded cylinder Cy = H x (0, (X) and study the effect of this truncation. We begin 
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with a result that shows the exponential decay of the solution to (12.81) ; compare with 
PSI Proposition 3.1], [531 Lemma 4.8] and |5SJ Proposition 4.1]. 

Lemma 3.1 (exponential decay). If Uq G tp £ ]HI®(11) D C(Cl) and f £ 

L^(0, T; L^(r2)), then, for every y >1, we have 

l|V'^llL2(o.T;L2(y<»_nx(r,oo)) 

where the hidden constant does not depend on neither ^ nor the problem data. 

Proof Consider, for a.e. t £ (0,T], the function w(t) £ that solves 

V-(y“Vr(;(t)) = 0 in C, vu{t)^g^c=0, trn w(t) = trnon 11 x {0}. (3.1) 

Since solves the fractional parabolic obstacle problem (12.81) and problem dUD has a 
unique solution, we immediately conclude that for a.e. t £ (0, T], we have w(t) = ‘f/{t). 
We now apply the decay estimate of [531 Proposition 4.1] to problem (EH) to obtain 

l|Vw(t)||L2(y._Ox(r.oo)) < e“^/^|]trn'^(<)||H'»(n)- 

Finally, integrating over time, and invoking the trace estimate (12.21) and the stability 
estimate (12.91) for problem (12.81) in terms of Uq, / and ip, we arrive at 

IIVr(;|]^2(Q 7’.^2('yQ f2x(y,oo)) e II^IIl2(q S), 

where 2) is defined in (1571) . This concludes the proof. 0 

The exponential decay of Lemma 13.11 allows us to consider a truncated version 
of the variational inequality (12.81) . To write this problem we define, for fk > 1, the 
Sobolev space 

= H\y^,Cr) : n;|r, = 0} , 

the convex set of admissible functions 

‘K.iCy) = ,Cy) : tTnw{x') > tp{x') a.e. x £ 

and the bilinear form 

ay{w,(j)) = ^ [ y°‘\7wV(j)dx' dy € Hy{y^,Cy). 

ds JCy 

With these definitions we consider the following truncated problem: Find v : [0,7] — 
“KiCy) such that tro ?;(0) = ug and for a.e. t £ (0, T) and all </> £ )^(Cy) 

(trn dtv{t),trn{v{t) - (j)))L 2 (^Q) + ay{v{t),v{t) - (p) < (f, trn(u(t) - ())))i 2 (n). (3.2) 

The analysis of this problem follows that of (12.5L developed in H2.3I For brevity, we 
only present the energy estimate 

II trn ^^lli~(o.r;L=(n)) + 11^11^2(0,T;#i(y°.c,r)) ^ 

We define Ha : IHI®(f2) Hf(y°‘,Cy), the truncated a-harmonic extension oper¬ 
ator as follows: if w £ in'*)!!), then W = HaW £ Hf{y°‘,Cy) solves 

V-(y“VW) = 0 in Cy, W = d on dyCy^jn x {‘j}, W = w; on 11 x {0}. (3.4) 
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We recall [551 Theorem 2.7] for problem (13.41) . If w G then 

II+ WdyyWh^y.^c,) < lkllHi+=(n). (3.5) 

The following result shows that by considering (13.21) instead of (12.81) we only incur 
in an exponentially small error; compare with [261 Lemma 3.3], [251 Lemma 4.3] and 
[24l Proposition 4.20]. 

Proposition 3.2 (exponential error estimate). Let and v be the solutions of 
(12.81) and (13.21) . respectively. Then, for 9 ^ > 1, we have 

II tl'a('^ — T)|||oo( 0 _T.i 2 (Qp + ||v('^ — 'c)||i 2 (o,T;L 2 (y“,C 5 .)) ^ ® 

where the hidden constant does not depend on neither , v, nor the problem data. 

Proof. By a trivial zero extension we realize that the solution v to problem (13.2|) 
belongs to then we can set ^ = u in (12.81) . We would like to set (j) = ^ in (13.2|) 

but, although it satisfies the constraints, it is not an admissible test function, as it 
does not have a vanishing trace at y = fT. For this reason, instead, we set cj) = in 
(15151) . where p G IF^’°“( 0 , cx)) is the following smooth cutoff function: 

p{y) = l 0 <y<|^, p[y) = ‘^ {<y - y) < y < 9 ^, p{y)=Q 7>y- 

With these choices of test functions we add the ensuing inequalities to obtain 

f dt||trn('^ -^^)lli 2 (a) + - v) < ay{v,p'^ -v). 

We now notice that 

a('^, -v)= - u) + II V'^||^ 2 (yc_Qx(r.oo)) > - v), 

so that we obtain 

id 4 ||trn('^ -v)\\l 2 (^Q) +ay{‘^ -v,'W -v) < ay{v, {p - 1)^). (3.6) 

It remains then to bound the right hand side of (13.61) . A straightforward computation 
reveals that if y < 9^/2 we have that {p — 1)^ = 0 , otherwise 

|V(y - < 2 + IV'^p^ , 

and thus 

II V(y - l)'^|li 2 (,.,c^) - ^ X /2 ^ C 2 L ■ 

Invoking a version of the Poincare inequality (12.11) based on the interval [ 9 ^/2, ‘j], we 
conclude IIV(p - < ||V'^||2,^^„_^^^^^2.y))- We now use this estimate 

in (13.61) and integrate in time to obtain 

II tra('^ - u)(t)||^2(Q) + f ||V(^ - u)(s)||^2(j^._cr) 

^0 

< II V^^llL=( 0 .T;L 2 (y“,Cy)) ll^'^llL2(0,T;L2(i/“,Ox(y/2,cx3))) ^6 

where the last inequality follows from Lemma [3Tl the fact that trn(^ — u)|(=o = 0 
and the stability estimate (13.31) of ([3.21) in terms of (D. Since t is arbitrary this implies 
the desired estimate. 0 
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4. Time discretization. We now proceed with the time discretization of prob¬ 
lem (HU. We could directly apply a suitable time discretization scheme to either (12.51) . 
(12.61) or (12.81) and then argue, for instance, by using the results of [55]. However, with 
the exponential convergence result of Proposition l3.2l at hand we will consider the time 
discretization of the truncated problem (13.21) based on the implicit Euler method. 

The discrete scheme computes the sequence C !fC{C^), an approximation to 
the solution to problem (13.21) at each time step. We initialize the scheme by setting 

trnT° = uo, (4.1) 

and for fc = 0 ,... /C — 1, let € “KiCy) be such that, for every (j) € X{Cy), 

ftro^^^,tro(H'=+i -</>)) + ay{V’^+\V’^+^ - cj)) 

\ /L2(f2) (4.2) 

<(f'“'+\tro(T'=+i-<(>))z,2(a), 

where 0 is defined in Section [5] and Itk^^ f dt G 

Existence and uniqueness of a solution to (14.21) follows directly from standard 
arguments on variational inequalities [gill]. The approximate solution to problem 
(HU is then defined by the sequence C IHI®(f2) where 

t7" = trnEE (4.3) 

Remark 4.1 (locality). The main advantage of scheme ()4.1I) - (I4.2|) is its local 
nature, which mimics that of problem (12.8|) . 

Let us now show the stability of the scheme. 

Proposition 4.2 (stability). Assume that uq G and f G L^(0,r;L^(H)), 

then € L°°{0^T] L^{n)) and V'^ € L^(0,T-, ,Cy)) uniformly inr. 

Proof. Set (j) = 'Ha'f’ in (14.21) . where T-Ca is the a-harmonic extension operator 
introduced in (13.41) . Upon denoting — Hafp we obtain 

V /L2(o) 

{f+\trn lU'=+i)i2(n) + ay{n^i^, 

The Cauchy Schwartz inequality and summation over k yields the result. □ 

The error analysis of (I4.1I) - (I4.2I) follows from the general theory presented in [28] . 
To present it we introduce the error 

E{v, V") = II trQ(?; - V'')\\l°^(o.T-,L^{Q)) + ||V(u - U'^)||L 2 ( 0 ,T;L 2 (y“,Cy))j (4.4) 

where V'^ and V'^ are defined in Section [2] We also define 

£{u,u^) = ||u — U’^||L“(0,T;L2(n)) + ||u — C’^||L2(0,T;H*(n)). (4.5) 

Corollary 4.3 (error estimates in time 1). If Uq G ^{£l) and^ G L^(0, T; L^(H)), 
then the solutions v of (|3.2I) and of (HU-(HU satisfy the uniform estimate 

^(^i ^ (l|uo||H'=(n) + ||f||L2(o_T;L2(n))) . 
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If, on the other hand, we have that Uq £ n]HI^®(r2), f £ BV{0,T; L'^{fl)) and 

= f_|_(tfc), then 

E{v, V^) < r(||f+(0) - (-A)«Uo|L.(o) + Vari=(n) f). 

In these estimates the hidden constants do not depend on v, nor the problem data. 
Proof. See Theorem 3.16 and Theorem 3.20 of [55]. □ 

The following result combines Corollary 14.31 the Caffarelli-Silvestre result [5], the 
trace estimate (12.21) and the exponential error estimate of Proposition 13.21 

Corollary 4.4 (error estimates in time II). Assume that Uq £ and f £ 

L^(0, T; then the solution u of (12.61) and the approximation U'^ defined by 

(ITSI) satisfy the uniform estimate 

£{u,U'^) < (||uo||H'»(n) + l|f||L2(o,T;L2(n))) + e 

If we have that Uq £ XiPi) D (11), f £ BV (0, T; L‘^{£1)) and f^ = f_|_(tfc), then 

f(u, Un < r(||f+(0) - (-A)®Uo||^.(^) + Vari.(o) f) + 

where in both estimates the hidden constants depend solely on the problem data. 

Proof. The definition of £(u, IP^), given in (14.5|) . the Caffarelli-Silvestre extension 
result trn '^ = u BM and estimate (12.21) yield £’(u,[/'^) < ,V'^), where E is 

defined in (|4.4D . Notice now that E is sublinear in its first argument, so that 


£{u, IP^) < E{v, V^) -I- II trn('^ — n)||L=x=(o,T;L 2 (n)) -I- ||^ — ?^||L 2 ( 0 ,T;L 2 (ycj_C 5 -))j 
The result now follows from combining Corollary 14.31 and Proposition 13.21 □ 

5. Space discretization. The results of previous sections are important in two 
aspects. First, we were able to replace the original parabolic fractional obstacle prob¬ 
lem m ( or any of its variants discussed in i i2.3l) by an equivalent one that involves a 
local operator (12.81) and is posed on the semi-infinite cylinder C. Then, we considered 
a truncated version (lO) of our problem, that is posed on the bounded domain Cy, 
while just incurring in an exponentially small error in the process. This is important 
because we shall discretize in space using first-degree tensor product Hnite elements. 
Section Ul presents a first order discretization in time and applies the general theory 
of discretizations of nonlinear evolution equations |28j to provide an error analysis. 

It remains then to discretize in space and to study its effect. We will follow 
[24l [25l [26] . where it is shown that U, solution of (12.3L possesses a singularity as 
?/ j. 0, so that the use of anisotropic meshes in the extended direction y is imperative 
if one wishes to obtain a quasi-optimal approximation error. The latter combines 
asymptotic properties of Bessel functions with polynomial interpolation theory on 
weighted Sobolev spaces m, which is valid for tensor product elements that exhibit 
a large aspect ratio in y. These references also show how to exploit the tensor product 
structure of Cy to design such a mesh. For convenience we recall this construction. 

Let ipQ = {K} be a conforming and shape regular triangulation of ft into cells 
K that are isoparametrically equivalent to either a simplex or a cube (TT] |T5] . We 
denote by an the shape regularity constant of iPn- Let ly = {1} he a, partition of 
[0,fX] with mesh points 





J = 0, 




7 > 


3 

1 — a 



Vj = 


7 


(5.1) 
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We then construct a mesh of the cylinder Cy by ®Xy, i.e., each cell T S 

is of the form T = K x I where K € IJq and I G Xy. Notice that, by construction, 
=ff£^y = When ^ is quasiuniform with ~ we have ss 

and, if = max{diam(iir) : K G ^}, then M 

Having constructed the mesh 3^y we define the finite element space 

V(^y) := [W G C'°(Cy) : W\t £ V{K) ® Pi(/) VT e STy, W|r^ = 0} , 

where, if K is isoparametrically equivalent to a simplex, V{K) = Pi(-?sr) i.e., the set 
of polynomials of degree at most one. If if is a cube 7^(if) = Qi(if), that is, the set 
of polynomials of degree at most one in each variable. 

We remark that, owing to (15.11) . the meshes £^y are not shape regular but satisfy: 
if Ti = Ki X Ii and T 2 = K 2 x I 2 are neighbors, then there is cr > 0 such that 

hi^<(7hi^, hi = \I\. 

While this is crucial to capture the singularities present in the solution, it also requires 
the development of anisotropic error estimates on Muckenhoupt weighted Sobolev 
spaces as detailed in [26l [27] . 

5.1. Fully discrete scheme. To describe and analyze the fully discrete scheme 
we must introduce an interpolation operator that preserves positivity of traces. In 
what follows we assume that there is an operator H^y : L^(Cy) —> V(J^) that verifies: 

• Locality. If w is such that tro w makes sense and v' is a vertex in , then 
n 5 yic(v', 0) depends only on the values of w in an H-neighborhood of v'. 

• Positivity preserving. If w is such that trn w makes sense then 

tro rc > 0 > 0. (5-2) 

• Stability. If w is a-harmonic and w G H'^{y^,Cy) D fl\{y^,Cy) then 

^ + lklbi(y.,c^) (5-3) 

• Approximation. If w is a-harmonic and w G H^{y^,Cy) D Hy{y°‘,Cy) then 

= (5-4) 

• Superapproximation. For 6 G [0, 2] we have that if tro w G then 

||tro(u;-n5yu;)||i2(f2) <||u;||H'5(n), (5.5) 

where the hidden constant is independent of <5, a-nd w. 

An example of such a construction is presented in Section jB] 

Let us now describe the scheme. We define 


£ V('^) : trn IF > trn Hjjy'Hat/’ a.e. 12} , 

where SLa is the a-harmonic extension operator introduced in (13.41) . The fully discrete 
scheme computes the sequence C t}(_{Siy)^ an approximation of the solution to 
(13.21) at each time step. We initialize the scheme by setting 

~ n^y'HaUo. 


(5.6) 
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For fc = 0, • • • , /C — 1, G solves 

/ 5^^+^ \ 

tra ' -W)\ + aAV^+\V^+^ - 1^) 

V / L 2 (n) 

< {f’^+\tTn{V^+^-W))^^^^^ VVF G (5.7) 

Standard results on variational inequalities yield existence and uniqueness of 
for fc=l,...,/C — 1. To obtain an approximate solution to the parabolic fractional 
obstacle problem (El), we define the sequence C ]HI®(n) by = tro 

Remark 5.1 (properties of the scheme). The main advantage of (I5.7D - (I5.6|) is 
that provides an approximate solution to the fractional obstacle problem ED based 
on solving the local evolution variational inequality (15.71) . Its implementation is simple 
and requires standard components of a hnite element algorithm. 

We note that, if Uq G ]HI^+®(n), the continuity (15.31) of the operator implies 

< l|uo||H^+^(n). (5.8) 

Indeed, the regularity results of [26l Theorem 2.7] show that, if u> G ]HI^+®(n), then 

dyyHaW G L'^{y^,Cy) and Wx''HaW G L?‘{y°',Cy)] see also (13.51) . 

Let us now present an error analysis for El- We will do so under different 
assumptions on the problem data, thus obtaining different rates according to the 
smoothness properties of the solution. We introduce the errors 

e- = e"(R-,F4) = I/--F4, (5.9) 

where R’’ and solve the problems (I4.1D - (I4.2I) and (I5.6I) - (I5.7D respectively, and 

=£^(V^) = v^ (5.10) 

5.2. Analysis with minimal regularity. Here we only assume that the right 
hand side satisfies f G L^(0, T; L^(H)), -0 G ]HI®(fl) and 0 < uq G ]HI^+'*(H). As a first 
result we obtain an a priori estimate for 

Lemma 5.2 (a priori estimates on V^^). If f G L^(0, T; L^(H)), 0 G ]HI®(H) and 
0 < Uq G ]HI^+®(H), then the sequence solution to (15.61) - (15.71) . satisfies 


II tra5V]^y|li2(Q) + T||VV)^||^2(y<»^Cy) + ll^^^.y^llL2(y“,Cy) 


fc=i 


fc=i 


5, T 


i2(o,T;L2(a)) + ||uo||Hi+'»(n) 


where the hidden constant does not depend on neither v, nor the problem data. 

Proof Set W = in (15.71) and multiply the obtained result by r. Using the 
identity 2a{a — 6) = — 6^ + (a — 6)^, with a, 5 G M, we derive 


tro t)U^+^|||2(Q) + — i>||VU^f ^|||,2(y<,_c,^) + ||Vf)U^f^|||,2(y,,^C;x) 


< 


r||f'=+'llL^(a)l|troi)U^+i|U2(n)<-r2||f 


■fe+l||2 


L^{n) + 9II trai>U^f^|li2(Q). 
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Adding this inequality over A: = 0,..., /C — 1 yields 


2 X] II + 2 II^K^||i2(yc C;x) 


2 X] ll^^^^rlli='(y“.Cr) 


fe=i 


< 


r 

2 


l|f|li2(0.T;L2(n)) + Iivr^jii2(,.,c,) 


The assumptions on f and Uq imply, in light of (I5.8|) . the asserted estimate. □ 

With these a priori estimates we can provide a first error analysis. Notice that by 
means of the change of variable u ■(— u — '0 one can assume that the obstacle 0 = 0. 
In this case then we have that C ac(Cy). 

Theorem 5.3 (error analysis with minimal regularity). //f £ L^(0, T; i^(0)), 
0 = 0 and 0 < Uq £ then 


II fro e’’|||cx=(o^T;L2(a)) + IVe^lli2(o,T;L2(y“,C5.)) ^ 

II fro e°|||2(n) + l|V^||j^2(o^T;L2(y“,Cy)) ’ 

where the hidden constant depends only on ||f||L 2 (o_T;L 2 (n)) ||uo||Hi+®(n)• 

Proof. Set 0 = in (14.21) and W = in (15.71) and add the resulting 

inequalities to arrive at 


(trot)e'=+\trne'=+^)^ 2 (o) + ,c.y) < 

- (trat)iyy,tra£'"+^') - Tay(iyy + r tro 

V J L^{Q) J'' \ J L^(Vt) 

where e^ andare defined by (15.91) and (I5.10L respectively. Added over k = 0,... 
1 , this inequality yields 


II tro e^llL2(n) + '^X ll^®^lli=(y“,Cy) ^ II ti'ne°y2(Q) + 


fc=i 


__i 
T 2 


' V 

X II lli^in) ] 

.k=0 ) 


1 

2 


+ \\^"^^y\\L'^{0,T-L‘^{y<^,Cr)) + l|f||L2(o,T;L2(n)) 


X 


IVe^||L2(o,T;L2(y<»,C:x))- 


Notice now that Lemma lOl implies that T~^J2k=o II ^Xy llL2(n) ^ ^ 
llVty \\L°°(o^T-L^{y'^,Cy)) ^ 1- These estimates allow us to conclude. □ 

Remark 5.4 (suboptimal estimate). Notice that, in the conclusion of Theo¬ 
rem [531 while the terms on the left hand side are squared, the interpolation error on 
the right is not. Therefore, even if the operator exhibited optimal approximation 
properties, this estimate is suboptimal in space. Nevertheless, from this result one 
can conclude convergence for rather general initial data uq and forcing term f . 

5.3. Analysis with regularity. The results of [7] show that, if Uq = 0 £ C^{^) 
with (—A)'*0 £ C'°’^“®(f2) and 0 < f £ T], C'°’^“®(fl)), then u satisfies: 

diU, (-A)^u £ logLip(L!)((0, T], C°’i-^(n)) s < i 


(5.11) 
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and 


diU,(-A)^ueC°’^((0,T],C°’i-*(ii)) s>i (5.12) 

These results, however, account only for the regularity of u, i.e., the regularity of 
tro ■ For the elliptic obstacle problem [51 Theorem 6.4] studies the regularity of the 
solution over the cylinder Cj and, on the basis of their findings, we shall assume that 

s < i ^ ^ e s> (5.13) 

Let us now present, under these improved regularity conditions, an error analysis. 

Theorem 5.5 (analysis with regularity). Assume that Uq, '4’ f are such that 
(15.111) (15.131) hold for ^ and uniformly in r. If ff^iy « then we have 


E{^,Vfr,)<r + \\ogMY 


M ^ + ||V£||Loo(o,r;L2(ya_c,^)) + 


ri/2 


where the hidden constant depends only on the problem data. 

Proof The results of Proposition 13.21 and Corollary 14.31 reduce the analysis to 
estimating the difference e'^ = . Using the well-known indentity 2a(a — b) = 

— b'^ + {a — by we derive 


- tro -I- II tro Ili 2 (n)) + '^W'^,ay) 

= (trn , trn ^ || Ve'=+i . 

Therefore, invoking (15.9p and (|5.10p . we arrive at 
(tro5e'=+\trne'=+i)^,(^^ = (tro i)e'=+\ tro 

-f (trnc)e'=+\trn(n5yU'=+i < ^|| tro |li2(n) 

+ ^I|trae''+^||i2(n) -f (trn0e'=+\tra(nj7yU''+i - Uj-f^)) 


L2(n) 


and 


-||Ve'"''‘^||^2(yo_Cy) — ay{,e' 


,efc+i) + ay(efc+i,n^,U'=+i -U|+i 


< i^l|Ve'=+i||i2(,.,c,) + ^l|V£'=+^||i2(,.,c,) + a^(e'=+\n^,U'^+i - U^^i). 


2d 

Consequently, 


(gfe+i 


2^11 *^f2e^'''^|||2(n) -I- ^||Ve^+^|||2(j^,»_Cy) 

+ ftrac)e'^+\tra(n^^U'=+i - V^+y) + Tay{e^+\ll^^V’^+^ 

\ ” / (f2) 
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Since £ ^(^), we use the scheme (15.71) with W = to derive 

‘£"+^ < l|troe'=+^||i2(n)+T||V£'=+^||2,(^„_^^^ + 

(tra5F'=+\tra(n^^F^'+i - V^+^)) + Ta^(V'^+\ns,V>^+^ - V^+^) 

The smoothness assumptions on allow us to integrate by parts to obtain 
ay(F'=+\n^^y'=+i-y|+i) = -l f Viy<^VV'^+^)in^^V'^+^ -V^+^)dx'dy 

S J C'y 

+ ^ (d,o.V^+\iTn^^,V^+^-trnV^+^) 

ds V ^ '^5' /L2(f2) 

Since = T“^tro0F^+^ — + d~^di,aV^^^ > 0, the regularity assumptions 

imply, in particular, that Z'^ £ (^([O, T]; C'°’^“'*(f2)). Therefore 

< l|tra£'=+i||i.(o) +^l|V£'=+^||i.(,..c,) 

+ T f trn(n^^F'=+i - V^+^) dx'. 

Jnx{o} 

We proceed now as in PH Theorem 4.24] and realize that it suffices to consider 
^ [ Z'=+i(trn(n^^F'=+i-n^^HaV')-tra(F'-'+i-7^„V’))dx'= ^ I{K). 

We analyze separately the cells K G according to the value of trn(17^+^ — 'Ha4’). 

• trn(l^^'''^ — 'Ha'tp) > 0 m a neighborhood of K. In this situation Z^+^ = 0, and 
thus I{K) vanishes. 

• trn(I^'''+^ - 'Ha'ip) = 0 in a neighborhood of K. The linearity of Il^y yields that 
I(A') vanishes. 

. tra{V’^+^ - is not identically zero nor strictly positive in a neighborhood 

of K. In this case, either ir^fUsTyV^^^ — Tl^y'Hafj) = 0 or Z^+^ =0. If i£ 
is such a cell, then there is Xq £ i£ where trnV^^^(xo) = '0(xo) so that the 
smoothness assumptions on ip and the regularity results of [7] allow us to conclude 
the growth estimate 0 < trnl^^^^(x') — '0(x') < Vx' £ K. By the same 

reasoning 0 < Z^+^(x') < Then, I{K) < hf^^. 

Collecting the derived estimates we obtain 

^ II trne^^^||i2(n) + fi^) +Th%-^. 

Add this expression over k. Using Proposition I4.2I and the regularity results (I5.1ip 

and (I5.12I) we have that tro U"^ £ L°°(0, T;]HI^+'*(U)). Therefore, using the superap¬ 
proximation of traces of the operator 11 (1531) we obtain 

K. I ^2(l+s) 

II trn £^'|li2(n) ^ “11 trn£^’'llio°(o,T;L2(a)) ^ £ ■ 

fc=i 

To conclude, we recall that h^r^ « □ 
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6. Positivity preserving interpolation over anisotropic meshes. The er¬ 
ror analysis that was presented in previous sections relied on energy arguments, and 
for that we needed to choose suitable test functions in the semidiscrete and discrete 
schemes. This brings forth the need for a positivity preserving interpolant. 

The construction of positivity preserving approximation operators has a rich his¬ 
tory in approximation theory. The classical convergence |19j and impossibility [20] 
results of P.P. Korovkin come immediately to mind in this respect. In the finite ele¬ 
ment literature, a positivity preserving interpolant was constructed in m and it was 
later showed in [50] that it cannot be of order higher than one. 

The operator of [10] is analyzed under the assumption that the mesh is shape regu¬ 
lar. In our setting we need an interpolant that preserves constraints of traces and that 
exhibits suitable approximation properties in weighted spaces and over anisotropic 
meshes. This makes the extension of the ideas of m not straightforward, if at all 
possible. For this reason, we will restrict our attention to the case s > 3/8 and, by 
combining the ideas developed in the construction of the positivity preserving inter¬ 
polation operator of [ 10 ] with the quasi-interpolation operator analyzed in [MIEZ], 
construct a positivity preserving operator on anisotropic meshes that possesses suit¬ 
able approximation properties on a-harmonic functions. 

Let us present a slight modification of the quasi-interpolation operator £ 5 -^ of 
[26l[27]. To do so, we introduce some notation and terminology. Given we denote 
by Af the set of its nodes and by A/in the set of its interior and Neumann nodes. For 
each vertex v € A/", we write v = (v', v"), where v' corresponds to a node of 5^, and 
v" corresponds to a node of the discretization in the extended dimension. We define 
hyi = mm{hK ■ K G v' 9 K}, and h^n = min{ft,/ : / € Xy, v" 9 /}. Given 
V G A/”, the star or patch around v is defined as = UtsvX, and, for T G we 
define its patch as St = UvgtS^. We set A/n = {v' : (v', v") G A/jn}. 

Let fii G be such that supp^i C where Br denotes the ball in 

R" centered at zero and with radius r < l/un; moreover, we require that p,i > 0 , 
f Ax' = 1 and that /ii has vanishing first order moments, i.e., / iii{x')x'' Ax' = 0 

for all f = l,...,d. Let ^2 S C'cr(®) be such that supp/r 2 C (0,ry), where < 
1/ct and / pi 2 {y)Ay = 1. We then define fi{x',y) := fii{x')fi 2 {y), which satisfies 
f y. Ax' Ay = 1 and supp y C Br x (0, r^)- For v G A/in we define 

, 1 f x' -v'\ 1 

y-iy{x ) = y J, y-2,v"{y) = j^y2 

and yv{x',y) = yi.y{x')y 2 y’{y)- We note that supp y^ C and fg y^ Ax'Ay = 1 
for any node v G A/in. 

Given a function w G L^(Cy) and a node v in A/in, the regularized Taylor polyno¬ 
mial of first degree of w about v is defined as follows: 

Wv{z) = / P{x,z)y^{x)Ax, 

Js, 

where P denotes the Taylor polynomial of degree one in the variable 2 : of the function 
w about the point x, i.e., P{x, z) = w(x) -I- 'Vw{x) ■ {z — x). 

If Av denotes the Lagrange basis function associated with the node v in the 
discretization SAy, we then define the averaged interpolant as follows: 

^ Wv(v)Av 
veA/in 



( 6 . 1 ) 
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is linear, stable and possesses optimal approximation properties in ,Cy). 

It is well suited for anisotropic meshes. The only difference between this construction 
and that of [13 [57] is in the particular choice of the weighting function fi. 

We now turn to the construction of m and slightly modify it to suit our purposes. 
Given cj) £ L^{n) we define the interpolation operator —>■ V(=%) by 


T^3rn'P= X] 

v'GA/o 


fj.iy{z')cj){z')dz' Av 


( 6 . 2 ) 


where 5', denotes the star (in fl) around v' and Av' denotes the Lagrange basis 
function associated with node v' in 5^. The assumption that fJ-i >0 yields the, fun¬ 
damental, positivity preserving property (15.21) . In addition, the symmetry properties 
(vanishing moments) of the function yiy imply that this operator preserves linears 
locally, therefore (15.51) holds. 

We now define a positivity preserving interpolant : Wi{Cy) V(5^). Let 
w £ Wl{Cy) and (v',v") be a node of then 

C 3 ryw{v',v") v'> 0, 

nyMv',0) v" = o. 

The approximation properties (15.41) of the operator 11are as follows. 

Theorem 6.1 (interpolation estimate). Let s > 3/8. If w G tly{y°‘,Cy) (~l 
H'^{y^,Cy) is a-harmonic and tiQW £ IHI^“''®(r2), then 


IlSyW{v' ,v") = 


||V(n; - IVyyW)\\L.^yyc,) < (#'^y)”^ (6.4) 

for 6 < 00 = min{l, The hidden constant blows up as 9 f 6o, is independent 

offf'^r and depends on w through \\wyy\\L2(^yfi^Cr)> \\'^^x'w\\L2(yc^Cy-), \\^M\L^iy<^,Cr) 
and II trn w||Hi+“(n)- 
Proof Consider 


||V(W - n 5 yw)||| 2 (yc_Cy) = ||V(n; - n 5 yw)||| 2 (j,a^Ox( 0 ,yi)) 

+ ll^(^ “ n5yU>)|||2(yoy2x(yi,r))> 

where yi is defined by (EB. Since over x (yi,lX) the operators and Il^y 
coincide, we invoke [13 Theorem 5.4] to arrive at 

||V(n; - Ilsyw)h2^yyaMyuy)) < (6.5) 

where the hidden constant depends on the function w only through its \\wyy\\L'^(yP 
and \\'\/\/x'w\\y 2 (^ya (;y'j norms. We then need to estimate the remaining term ||V(r(; — 
^^xW)\\L2{y‘^,nx(o,yi))- To do this, we proceed as follows: 

l|V(w-n5yW;)||i2(j^.^Ox(0,yi)) < l|V(w-£5yw)||L2(j,a_nx(0.yi)) 

+ l|V(£j?yW - n5yU;)||L2(y<.^Qx(0,yi)), 

where £ 5 ^ is defined as in ()6.1I) . The first term of the expression above is controlled 
by the right hand side of (16.51) by invoking, again, [13 Theorem 5.4]. To estimate the 
term — II over the first layer LI x (0, ?/i), we use the definitions of £ 5 ^ and 
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given by (16.11) and (16.3L respectively and exploit the fact that, for every node 
v' e we have Cs^w(v',yi) = ,yi), to write 






KeS^n ^ ^ v'GAfn 


where TZsn is defined by (16.21) . Av' denotes the basis function associated with the node 
v' in the discretization IJq and Ao(t/) the one that corresponds to the node t/o = 0 in 
the discretization ly defined by the mesh points (15.11) . 

Using the well known finite intersection property of the supports of the basis 
functions Av' we see that, to conclude, it suffices to estimate, for i = 1 ,..., d + 1, 


£i-=[ y°‘ [ \{C 3 nyw{v', 0 )-nsnw{v',0)) d^i{Ay,{x')Ao{y))\^ dx'dy. 


Jo J K 

To achieve this we notice, first of all, that if i = 1,..., d we have 



( 6 . 6 ) 


while, if i = d + 1, 



(6.7) 


\Jo / 

We must uniformly bound the difference between and TZ^^w over A/n- 

Since / y, 2 iy) dj/ = 1 and y,i has vanishing moments, we have that, for a.e. z' € 11, 


= f 

Js 


y-iy (x')w{x', 0) dec' 


L 


y-y,o){x',y) [w(x',0) + Va;''u;( 2 ', 2 /)(v' - x')] dx' dy. 


Using this we see that, for any v' G Afn, we have 



We now proceed to bound each one of these terms separately. 

• Bound on III: Using the scaling properties of the function yy,o) and the Cauchy- 
Schwarz inequality we obtain 
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Since w is a-harmonic and eu dictates the grading of the mesh on the extended 
dimension, we have that \\dyw\\i, 2 i^yc, where « #>%; see 

§5.2]. This yields 


|III| 


( \ ^1 /2 


Bound on I: A mean value result allows us to write 


|I| < / fJ.^^>^o)\dvwix',r]{y))\ydx'dy. 

‘^(v',0) 


After this, we proceed as in the bound for III to conclude 

1/2 


< 


1/2 


M-h (6.9) 


• Bound on II: We use that w is a-harmonic and a mean value theorem to get 
|n| < [ \Vl,wi(ix',z),y)\dx'dy 

<7 t4 — (/ y~"dx'dy\ ||v 2 ,u;||i 2 (j ^„,5 ). 

K' yi / 

Using the local regularity results of [551 Theorem 2.9], we derive 

ll^x''*^lli2(y°,nx(o,yi)) ^ yill 
This allows us to obtain 


nil < 


d-l 1/2 ij 2/ “dx'dy"! < 1/2 (/ 22 “dy^ . (6.10) 

^v' Vi V''‘®(v',o) / 2/i '^•'0 / 


To conclude we insert (|6.8D - (I6.10I) into (|6.6I) and (|6.7I) . If / = 1,..., d we have 
fz < y“ dy^ [I^ + Il2 + IIl2] 

ry^ \ r M-2 1 1 

= £i,l + ^z,2. 


< 




hi htr^ 


yi 


Using the fact that y“ G A 2 (M‘^“'"^) 

where we have used that M ?s rs (^/5^)1/('^+i) together with the fact that 

yi/h^' is uniformly bounded. We now proceed to estimate the term £i 2 - 

£^,2=(- r\“dy) r\-^dy) hiy,. 

\yi Jo ) \yi Jo ) 
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Since t/“ € and k, M ^ we have £i ^2 ^ M These yield £i < M 

Finally, we focus in the case i = d + 1: 

^d+i <^( r 2/“ dy) [f + If + III^] < M-2 + y-^ht,, 

Ui \Jo / 

where we have used, again, that y°‘ £ Since the grading parameter 7 

satisfies 7 > 3/(2s), we have that, for s > 3/8, do = min{l, ( 8 s — 3)/(4s)} > 0. 
Therefore, we obtain that, for 0 < 6 q, we can bound £d+i ^ 

Collecting the derived estimates for Si, with i = 1, • ■ • , d and Sd+i gives us (EH). 0 
We now show the stability of the operator 11 5 ^ (15.31) . 

Corollary 6.2 (stability). In the setting of Theorem \6.1[ we have 

^ \\'^vv\\L'^{yf^ ,Cy) + l|VVx'U’||L2(y<._c,x) + II tl'OU’llHi+^n), 

where the hidden constant does not depend on SSy. 

Proof. The result follows as a simple application of Theorem 16. II In fact. 


,Cy) < l|V(w - n^yW)||i2(yc,_Cy) + \Ww\\l^{v°‘M 

and (16.41) yield the desired estimate. □ 

To conclude let us use this operator to obtain error estimates. 
Corollary 6.3 (error estimate for Assume that s > 3/8 and set 


do = -1, s £ 



8s — 




4s 


s e 



In the setting of Theorem \5.5[ if t 


{ffSScr)'^y then 


<|log#^rr(#^r)"^ (6.11) 

for d < do. The hidden constant blows wp as d f do, is independent of and 

depends only on the problem data. 

Proof. The choice of r and, depending on the value of s, a comparison of the terms 
(^^)“®/(‘^+i) and [ff.Sry)-A+A/U+i).j--i /2 right-hand side of the estimate of 

Theorem 15.51 yields the result. □ 

Corollary 6.4 (error estimate for u). Assume that s > 3/8. In the setting of 
Theorem \5.5\ and Corollary 1 6. dl we have 

f(u,C/5.J 


for d < do. The hidden constant blows up as 9 f 60 , is independent of ffSSy and 
depends only on the problem data. 

Proof. The desired estimate follows from Corollary 14.41 and the estimate (16.111) . □ 
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